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A change in boundary conditions (BC) from uniform Dirichlet to non-identical BC on the edges 
of a triangular billiard often brings about a dramatic change in quantum spectral fluctuations. We 
provide a theory for this based on periodic orbits and show that non-identical BC on adjacent edges 
can lead to a quantum splitting (QS) of periodic orbit families which results in a significant change 
in the form factor. Thus, the classical spectrum alone cannot determine quantum correlations. 
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Recent studies jl|,[| on fluctuation measures in chaotic systems seem to indicate that the quantum correlations are 
fully determined by the classical spectrum of the Perron- Frobenius operator ||. Using different approaches, Agam, 
Altshuler and Andreev Q| and Bogomolny and Keating g] show that the diagonal and off-diagonal part of the density- 
density correlation are related to a purely classical quantity, which under some approximation reduces to the classical 
t— I . zeta function j| . 

There are several fallouts of the AAA-BK theory. One that has been scrutinized recently by Prange |4| concerns 
possible deviations from random matrix theory results and the conditions under which this can be observed. Another 
consequence (and one that is of relevance here) is the absence of parity effects in fluctuation measures. In other 
words, the AAA-BK theory predicts that quantum systems having the same classical dynamics exhibit identical 
quantum correlations. There is however a tacit assumption in the BK approach which leads us to this conclusion - 
that degenerate periodic orbits with identical actions and stabilities also have the same quantum phase. There can be 
Q\ , instances however when this is not true. For example, arithmetic billiards abound in degenerate periodic orbits and 
exhibit Poisson fluctuations when the boundary conditions are Dirichlet j5). However, when the boundary conditions 
are not identically Dirichlet, pairs of degenerate periodic orbits can have phases differing by tt leading to a net decrease 
"-^ ■ in the form factor H . 

The effect of boundary conditions can be even more significant in planar triangular billiards and we shall deal 
with these henceforth. Of these, the ones that are integrable have internal angles of the form ir/rii and in all these 
cases, degeneracies exist in the classical periodic orbit actions of topologically distinct orbits leading to non-universal 
spectral fluctuations [7j||. Thus, the spectral rigidity As(L), increases with a slope larger than yg- JhJ in the 
region L << L max where L max is determined by the frequency of slowest oscillation in the quantum density p(E) 
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Generic rational triangles on the other hand have internal angles of the form irmi/rii where TLf =1 rrii ^ 1. These are 
referred to as pseudo- integrable (PI) billiards 12 iljj. As in integrable systems, their invariant surface in phase space 
is 2-dimensional but the topology is that of a sphere with multiple holes and not a torus . This difference leads to 
a rather dramatic change in the quantum eigenstates. The eigenfunctions for example often exhibit irregular nodal 
patterns and a Gaussian amplitude distribution p6[ while fluctuation measures display a behaviour [p]p^] that ranges 
from the integrable [Q to the chaotic |]l8|-^o|,[7| limits. 

There are several interesting questions concerning quantum fluctuations that polygonal billiards throw up. A point 
that has often been debated is the role of diffractive periodic orbits in determining spectral measures [^l| . Admittedly, 
the quantum spectrum does know about these orbits [ p2| though its importance in determining spectral measures 
is possibly negligible |23|| . A related question concerns the effect of boundary conditions on spectral fluctuations. 
To illustrate this, we refer to fig. 1 where the rigidity, A3 is plotted as a functions of L for the right triangle 
(7r/2, 7t/3, 7t/6) with (a) Dirichlet boundary condition on all the three edges and (b) Neumann boundary condition 
on the edges enclosing the right angle and Dirichlet BC on the third. One might argue that the crossover is related 



o 



1 

o 



journal ref: Phys. Rev. E, VOLUME 57, PAGE R3699, 1998 
email: dbiswas@apsara.barc.ernet.in 



1 



to the fact that case (a) is integrable while (b) is not |24|. We have thus verified that there is indeed a shift with 
boundary conditions in genuine pseudo-integrable enclosures such as the (tt/2, 37r/10, 7r/5) triangle. 

At first glance, it may seem that apart from an overall phase factor, the contribution of each periodic orbit family 
is identical for the two cases of the (tt/2, 7r/3, 7r/6) triangle. We shall demonstrate here that this is not the case. 
To this end, consider the family of periodic orbits shown in fig. 2. Orbits Fl and F2, belong to the same family 
classically though the quantum phase accumulated by these differ by n when edges 1 and 3 have Neumann boundary 
conditions (we refer to this as case (b) while (a) denotes Dirichlet BC on all edges). In other words, the family of 
periodic orbits split up in the quantum-mechanical sense in case (b) while case (a) preserves the full classical family. 
The semiclassical density of states, 
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p{E) ~ p av (E) + X / — J2Y1 ~W=T= cos(krl p - J - rn p ir) (1) 



for the two cases are thus distinct. In the above, E — k 2 while a p refers to the area occupied by a primitive periodic 
orbit family characterized by their length l p and the number, n p of reflections from Dirichlet edges. Eq. ([!]) has 
only the leading order fluctuation in the density of states. It neglects the contribution of isolated periodic orbits and 
diffractive periodic orbits. Also note that orbits that occur in families necessarily undergo even number of reflections 
from the edges so that the net phase in case (a) is zero. 

With this background, we are now ready to explore the effect of BC on spectral measures. For two point correlations 
such as Aa(L) and T,2(L), the central object is the form factor 4>(T) = f_ oo R2(x)exp(ixT/h)dx where R2{x) is the 
2-point spectral correlation function (R,2(x) — (p(E + x)p(E))). Expressed in terms of periodic orbits, 4>{T) — 
J2i J2j A i A 3 cos(5i - Sj)S(T - (Ti + Tj)/2)\ where A t = Cai/^/Ml for marginally unstable billiards, Si = kl t - 




tt/4 - n.TT, Ti = dSJdE and C = ^1/(32tt 3 ). 

It is customary to analyze the diagonal and off-diagonal part of 4>{T) separately and we first show that for case (b), 
the diagonal contribution 4>d(T) — Af5(T — Tj)) is smaller as compared to case (a). Let us assume, that the 
family labeled by i splits up quantum mechanically in case (b) into two parts occupying areas an and ai2 respectively 
where an + = QSj. Its contribution to 4>d(T) is thus proportional to a 2 x + a 2 2 while in case (a), it is proportional 
to a 2 x + af 2 + 2ai\Oi2- Further, since the two parts of the classical family have a different phase in case (b), there is 
an off-diagonal (OD) contribution from within this classical family. Its magnitude is proportional to 2anai2 cos(-7r) so 
that the net decrease in contribution of a single classical family is proportional to Aanai2- 

Note that the off-diagonal part of the form factor has cross contributions as well where parts of two distinct classical 
families are involved. When the classical dynamics is integrable and no QS occurs, the OD contributions average 
to zero in the absence of degeneracies amongst periodic orbit actions. Thus the diagonal contribution equals the 
asymptotic value of <j>(T) which equals p av /2-K. This asymptotic law is referred to as the semiclassical sum rule |fll| . 
Even in the presence of QS, the semiclassical sum rule holds. Thus, cross terms involving parts of distinct classical 
families do contribute in case (b). In summary then, the following comparison between cases (a) and (b) can be made 
when there are no degeneracies in the lengths of topologically distinct periodic orbits. For T << Tjj, the form factor 
equals 
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<t>(T) = {C 2 Y, TT^ T -Ti))..- case (a) (2) 
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=f 2 E |ai(2 ;: 1)}2 ^-^))--w (3) 

while in both cases, 4>(T) = p av (E) / (2w) as T — > oo. Here an — oiidi, = (1 — Q-ijcn and Th is the Heisenberg time. 

Note that in an integrable enclosure, the areas are identical for almost all orbit families so that for case (a), 
a straightforward application of the proliferation law for periodic orbit families leads to the conclusion that </>(T) is 
constant and equals p av (E)/2ir. For case (b) however, the factor (2a^ — l) 2 varies with the orbit and depending on 
the splitting mechanism, </>(T) may be explicitly T dependent even in an "integrable" enclosure ]25| ], 

In order to concretize these notions, let us take another look at the (ft/2, ir/3, tt/6) enclosure of fig. 1. For this 
integrable billiard, the length spectrum can be expressed in terms of winding numbers on tori and it is easy to verify 
that there exists degeneracies in the lengths of topologically distinct periodic orbits. For case (a) then, the sum in 
Eq.(Q) is over distinct lengths U instead of topologically distinct orbits. Correspondingly, the area <Zj should now be 
interpreted as the total area occupied by all orbit families having length Zj. An immediate consequence is that 4>(T) 
is no longer a constant for all T since the degree of degeneracy varies with length Q. A plot of J(r) = -y- J 4>(t')cIt' 
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with respect to r = T/(2np av ) is provided in fig. 3. For generic integrable systems without degeneracies in periodic 
orbit lengths, I(t) = t while in the present case, one observes a non-linear increase. 

For case (b), the splitting mechanism needs to be incorporated and for this example, the ratio in which certain orbit 
families split up has been arrived at by Shudo ^(|. As before, the sum in Eq.(Q) now refers to distinct lengths while 
the effective area ai(2cti — 1) (denoted by Hi) is the sum of all areas occupied by degenerate orbit families weighted 
appropriately by the phases. Thus 5j = ^2 k (— l) nk ak where ak is the area occupied by a family having length Zj and 
which undergoes reflections from the Dirichlet edge. 

Once more, rather than the asymptotic proliferation rate of periodic orbit families, it is the variation of Hi with 
length [^7j which determines the form factor. A plot of I(t) for case (b) (see fig. 3) reveals a non-linear increase 
having a smaller overall slope and a form that is distinct from case (a). Thus a change in BC from uniform Dirichlet 
to non-uniform BC leads to a significant change in the form factor. This difference indeed shows up in the spectral 
rigidity, A^(L). In fig. 4, we compare the predictions of periodic orbit theory with the exact (numerical) values of A3 
in the range 4 < L < 10 [^8| for cases (a) and (b). While the agreement for case (a) is excellent, the predictions of 
periodic orbit theory capture the overall behaviour in case (b). 

The discussion so far holds for all polygonal billiards where adjacent edges enclosing an angle of the form ir/rii 
have non-identical (NI) boundary conditions. In such cases, periodic orbit families do not split up at this angle 
classically though as demonstrated earlier, they can split up semiclassically. For angles of the form miir/rii (mi > 1) 
however, orbit families do split up classically and traverse different paths thereby reducing the extent of periodic orbit 
families. Thus different sets of boundary conditions only result in an overall phase factor for each family and hence 
do not significantly affect spectral measures. In exceptional cases however, the effect of boundary conditions can be 
significant. This can be observed when the angle is of the form jn.i (m^ > 1) but close to an integrable wedge. 
As an example, consider the rigidity A^(L) for the irrational triangle (n/2, tt /V9.1) which is close to the integrable 
(7r/2, 7r/3) enclosure (fig. 5). Case (a) clearly exhibits fluctuations close to Poisson while case (b) shows typical GOE 
fluctuations for the energy range considered. Note that the triangle has infinite genus though over short time scales 
(less than the Heisenberg time) the dynamics hovers around its integrable counterpart while even after 10 9 reflections 
from the boundary, parts of the constant energy surface remain unexplored. The two non-integrable acute angles 
however serve to split up periodic orbit families though the lengths of the resulting families remain close to that of 
the original family in the integrable enclosure. This subtle re-organization of periodic orbit families leads to Poisson 
fluctuations in case (a) since the degeneracies in orbit actions which exist for the (tt/2, tt/3) triangle get lifted in the 
(it/2, 7t/\/9.1) enclosure. On the other hand, when the lifting of degeneracies is accompanied by a difference in phase 
between two split families (case (b)), the change is significant and leads to GOE like fluctuations for the energy range 
considered | f29fl . 




In summary, we have demonstrated that a change from uniform Dirichlet to non-identical boundary conditions on 
the edges of triangular billiard can lead to significant changes in fluctuation measures. This can be observed when 
the angle enclosed by the edges with NI boundary conditions is of the form 7r/n or sufficiently close to it J30|. The 
mechanism involved is quantum splitting due to which adjacent families having (almost) identical lengths acquire 
different phases leading to a significant drop in contribution from both the diagonal and off-diagonal terms in the 
form factor. In particular, we have shown that it is possible to explain the spectral fluctuations of the (tt/2, 7r/3, 7r/6) 
triangle when the boundary conditions are not identically Dirichlet using periodic orbit theory. We conclude by noting 
that there exist quantum systems whose density correlations cannot be determined fully by the classical spectrum. 



[1] O.Agam, B.L.Altshuler and A.V.Andreev, Phys. Rev. Lett., 75, 4389 (1995); see also A.V.Andreev and B.L.Altshuler, 
Phys. Rev. Lett., 75, 902 (1995); A.V.Andreev, O.Agam, B.D.Simons and B.L.Altshuler, Phys. Rev. Lett., 76, 3947 
(1996). 

[2] E.B.Bogomolny and J.Keating, Phys. Rev. Lett. 77, 1472 (1996). 
[3] P.Cvitanovic and B.Eckhardt, J.Phys. A24, L237 (1991). 
[4] R.Prange, Phys. Rev. Lett., 78, 2280 (1997). 

[5] E.B.Bogomolny, B.Georgeot, M.-J.Giannoni and C.Schmit, Phys. Rev. Lett. 69, 1477 (1992); J.Bolte, G.Steil and F.Steiner, 

Phys. Rev. Lett. 69, 2188 (1992). 
[6] R.Aurich, F. Schemer and F.Steiner, Phys. Rev. E51, 4173 (1995). 
[7] A.G.Mittenburg and Th.W.Ruijgrok, Physica A210, 476 (1994). 
[8] D.Biswas, M.Azam and S.V.Lawande, J. Phys. A24, 1825 (1991). 
[9] M.L.Mehta, Random Matrices, Academic Press, N.Y., 1967. 




3 



[10] In generic integrable systems, As(L) = L/15 while for the Gaussian Orthogonal Ensemble (GOE) of matrices, A3 

increases logarithmically. 
[11] M.V.Berry, Proc. Roy. Soc. Lond. Ser. A400, 229 (1985). 
[12] P.J.Richens and M.V.Berry, Physics D2.495 (1981). 

[13] Irrational triangles are possibly ergodic [1141 . The ones accessible to computations are rational and hence pseudo-integrable. 
[14] A. Hobson, J. Math. Phys. 16, 2210 (1976); R.Artuso, G.Casati and I.Guarneri, Phys. Rev. E55, 6384 (1997). 
[15] The invariant surface is labeled by the genus, g where 2g is the number of cuts required to make the surface singly 
connected. 

[16] D.Biswas and S.R.Jain, Phys. Rev. A42, 3170 (1990). 
[17] A.Shudo and Y.Shimizu, Phys. Rev. E47, 54 (1993). 

[18] M.V.Berry and M.Wilkinson, Proc. Roy. Soc. Lond. Ser. A392, 15 (1984). 
[19] T.Cheon and T.D.Cohen, Phys. Rev. Lett. 62, 2769 (1989). 
[20] D.Biswas, Phys. Lett. A159, 28 (1991). 

[21] Y.Shimizu and A.Shudo, Chaos, Solitons and Fractals, 5, 1337 (1995). 
[22] N.Pavloff and C.Schmit, Phys. Rev. Lett., 75, 61 (1995). 
[23] D.Biswas, Phys. Rev. E54, R1044 (1996). 

[24] Case (b) corresponds to the even parity states of the tt/3 rhombus, a pseudo-integrable billiard with g = 2. 

[25] In generic pseudo-integrable enclosures, {at} are orbit dependent p3|| . 

[26] A.Shudo, Phys. Rev. A 46, 809 (1992). 

[27] D.Biswas, Pramana - J. Phys. 48, 487 (1997). 

[28] Due to limited number of orbits, the fit is not expected to be good for small L. 

[29] Though we observe GOE-like fluctuations in several cases, these quantum systems do not necessarily belong to any 
universality class. There is a strong energy dependence that manifests clearly in A3 but is not so significant in the nearest 
neighbour distribution. 

[30] There are exceptional cases when the number of reflections from Neumann edges is even for every family so that no quantum 
splitting occurs. An example is the (7r/2, 7r/3, 7t/6) enclosure with Dirichlet BC on edges 1 and 2 and Neumann on edge 3. 
In this case the spectral measures are close to case (a). 



4 



FIGURE CAPTIONS 



Fig. 1. As(L) for the (w/2,w/3,ir/6) triangle with (a) Dirichlet boundary conditions on all edges (A) (b) Neumann 
boundary conditions on the edges enclosing the right angle and Dirichlet on the third (o). The averaging interval is 
[e„ — Ae, e„ + Ae] with e„ = 800 and Ae = 300. The straight line and the smooth curve are respectively the Poisson 
(L/15) and GOE results. 

Fig. 2. Two periodic orbits Fl (thick line) and F2 belonging to the same classical family. 

Fig. 3. A plot of I{t) for the (7r/2, it/3, tt/6) billiard. The curve marked (integrable) is for case (a) while (PI) represents 
case (b). Also shown are three lines with slopes 1.0, 0.85 and 0.45 marked (c), (d) and (e) respectively. For averaging, 
see fig. 1. 

Fig. 4. The chain and dotted curves are the exact values of A 3 (L) for case (a) and (b) respectively. The diamonds 
and squares are estimates obtained using periodic orbits for case (a) and (b) respectively. For averaging, see fig. 1. 



Fig. 5. The rigidity for the irrational triangle (n/2, 7r/v9T). Case (a) exhibits Poisson fluctuations (o) while case (b) 
shows GOE fluctuations (+) when e n = 500 and Ae = 150. 
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